From recollisions to the knee: A road map for double ionization in intense laser fields 
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We examine the nature and statistical properties of electron-electron collisions in the recollision 
process in a strong laser field. The separation of the double ionization yield into sequential and 
nonsequential components leads to a bell-shaped curve for the nonsequential probability and a 
monotonically rising one for the sequential process. We identify key features of the nonsequential 
process and connect our findings in a simplified model which reproduces the knee shape for the 
probability of double ionization with laser intensity and associated trends. 
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Multiple ionization by intense laser pulses often takes 
place as an uncorrelated sequence of single ionizations- 
often, that is, but not always: In fact, it has furnished one 
most striking surprises in physics of recent yearshj when 
early experiments on double ionization [1, 0, 0, clearly 
showed the unexpected existence - and even dominance- 
of another channel -the so-called non-sequential double 
ionization (NSDI) channel-which cannot be regarded as 
two uncorrelated single- ionization events. What was even 
more surprising was the prominence of this channel: The 
by-now famous "knee" shape indicates that at some in- 
tensities, the yield from this process exceeds the contri- 
bution from the sequential channel by many orders of 
magnitude 0, d, 0, H[ . Knees have been observed in the 
yields of multiply charged ions too, suggesting the 
significance of nonsequential channels there as well. 

Today NSDI is regarded as one of the most dramatic 
manifestations of electron-electron correlation in nature. 
Despite two decades of intensive research, it precise quan- 
tal mechanism 0, HI, HI [ToL fll| remains an intriguing sub- 
ject. However, the gen eral outline of the process, known 
as "recollision" [Ta Il3j|. are clear enough [lij ]: The laser 
ionizes an electron, which picks up energy from the laser 
field and, upon reversal of the field, is hurled back at the 
core where it dislodges a second electron [l^, EH ■ What 
is less clear is the nature of these collisions: How efficient 
are they, how much energy is shared during the collision, 
and when do they lead to double ionization? After all, 
not every recollision leads to double ionization, or does 
so right away. These questions have a direct bearing on 
the NSDI probability, and of course, the resulting knee 
shape. 

In this Letter, we examine the collision process in de- 
tail and answer these questions using classical mechan- 
ics [HI, [H, [I?], El, EH , which is increasingly widely rec- 
ognized as the tool of choice for ionization phenomena 
sufficiently high above threshold, since, remarkably, clas- 
sical and quantum-mechanical collision results are virtu- 
ally indistinguishable in that regime This a gree ment 
is ascribed to the dominant role of correlation [lq |. Af- 



ter identifying characteristic features needed for double 
ionization, we separate the knee into two contributions 
which peak at differing intensities: A bell-shaped curve 
for the nonsequential process and a monotonically rising 
curve for the sequential component. As an ultimate test 
of our investigation, we connect the bare essentials of the 
collision process through a mapping which transforms ini- 
tial conditions to ionization probabilities. This mapping 
yields the hallmark bell shape for NSDI. However, this 
does not exhaust the information contained in the map- 
ping: Remarkably-and to quote just one example- the 
fall-off of the NSDI probability with intensity turns out 
to be not solely due to the depletion of the initial sample 
by sequential double ionization (SDI) but also due to the 
decreasing efficacy of electron-electron collisions with in- 
creasing collision energy, a fact that comes out directly 
from our collisional model, as we will show below. 

In order to investigate these issues, we consider the 
following Hamiltonian system [l|| describing a one- 
dimensional Helium atom using soft Coulomb potentials 
driven by a linearly polarized laser field of amplitude Eq 
and frequency u> : 
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where x and y denote the position of each electron, and p x 
and p y their (canonically) conjugate momenta. 

Without the laser field (Eq = 0), typical trajecto- 
ries associated with Hamiltonian ([1]) are composed of 
one electron close to the nucleus (the "inner" electron) 
and one electron further away (the "outer" electron) . 
When the laser is turned on, the outer electron quickly 
ionizes while the inner one experiences a competition 
between the laser excitation and the Coulomb interac- 
tion with the nucleus. Its effective Hamiltonian is given 
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by [191: 



p 2 2 

H in (y,Pv,t) = +yE smwt. (2) 

1 Vv + 1 

Typical trajectories of Hamiltonian H m display two dis- 
tinct behaviors : Close to the nucleus, the Coulomb inter- 
action is strong enough to keep the electron attached to 
the core (therefore referred to the "bound" region in what 
follows). In contrast, further away from the nucleus ("un- 
bound" region), the laser field prevails and sweeps the 
electron away. The separation between the two regions 
is very sharp which makes it easy to compute numeri- 
cally the size of the bound region as a function of the 
laser field intensity [l9(. As a result, at the beginning of 
the laser excitation, if both electrons are in the unbound 
region, one can reasonably expect sequential double ion- 
ization, while if at least one electron is in the bound re- 
gion, a recollision (with the outer electron) is needed to 
raise it to the unbound region in order to ionize. When 
both electrons ionize at (about) the same time, the pro- 
cess is usually labeled as nonsequential double ionization 
(NSDI), whereas double ionizations with a large delay be- 
tween ejections are sequential double ionizations (SDI). 
However, these definitions fail to take into account the 
correlated nature of multiple ionization processes. For 
instance, a recollision may put the inner electron into an 
almost-bound state which then takes a long time (some- 
times more than one laser cycle) to ionize. With the 
previous definition, these so-called "recollision excitation 
with subsequent ionization" (RESI) [g, [20j] events - by 
no means rare - would be labeled as SDI whereas they 
clearly correspond to a correlated process in the same 
way as NSDI. Rather than decomposing the double ion- 
ization yields into SDI and NSDI contributions, here we 
consider the decomposition into correlated double ion- 
ization (CDI), where at least one recollision is needed 
for double ionization, and uncorrelated double ionization 
(UDI) where no recollision is needed. For a good approx- 
imation of the UDI contribution to the double ionization 
yield we use the distance between the inner electron to 
the nucleus, obtained from Hamiltonian |(3J|. This dis- 
tance is best expressed in terms of its energy 
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The smaller this energy is, the closer to the nucleus 
the electron is. Since Hq is integrable, there exists 
(at least locally) a canonical transformation which maps 
this Hamiltonian into action and angle variables so that 
it only depends on the actions. The action of the 
inner electron is defined by A — <f> p y dy / (2ir) . We 
denote by A m (Eq) the action of the outermost in- 
variant torus of Hamiltonian |(3J). A good approxi- 
mation of the unbound region is given by T> (Eq) — 
{(y,p v ) s.t. A(y,p y ) > A m (E Q )}. The UDI probability 
is given by the proportion of the ground state energy 
surface where both electrons belong to T> (Eq). 
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FIG. 1: Double ionization probability for Hamiltonian (ft) (big 
black circles) and the mapping ([5]) (small red circles) as a 
function of the laser intensity I for a wavelength of 780 nm. 
We also indicate the expected UDI probability obtained as 
described in the text (continuous blue curve) . Vertical dashed 
curves refer to intensities where phase portraits are displayed 
in Fig. El 



In Fig. [U we display the UDI component predicted 
by this model (continuous blue curve) which is in very 
good agreement with the double ionization probability 
obtained by integrating the full Hamiltonian |T]) in the 
high-intensity regime. 

Statistical analysis of recollisions. From the previous 
definition of UDI, the CDI probability is obtained by sub- 
tracting the expected UDI from the double ionization 
yield. 

A quick inspection of Fig. [T] reveals a bell-shaped curve 
for the resulting CDI component. A rather intuitive 
mechanism to explain the decreasing part of this bell 
shape is a conversion from CDI trajectories into UDI, 
when the laser field becomes stronger. However, in Fig.[TJ 
we notice a local decrease of the total yield (also observed 
with quantal computations [2l[) which is larger by several 
orders than the increase of UDI. We notice that this in- 
compatibility is readily observed in Fig. 1 of Ref. [2l|]. As 
we shall see below, the decrease of CDI is mainly due to 
the decrease of recollision efficiency with the laser inten- 
sity. In order to identify the key features responsible for 
the decrease of CDI, and thus for the knee shape, we first 
collect statistical data from the trajectories associated 
with Hamiltonian ([I]) at recollision times. In particular, 
we record the number of recollisions, the momentum of 
the outer electron at recollision and the energy (or rather 
the action) exchanged during non- ionizing recollisions. 

It is well-known that the maximum energy the outer 
electron can bring back to the core is equal to [l^, HH 
£max = kEq/ (4a; 2 ) where n « 3.17. However, an inspec- 
tion of the mean energy exchanged during recollisions 
shows that this amount is significantly smaller and de- 
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creases with intensity during (non-ionizing) recollisions 
(see Fig. ^ at relatively high intensities. This is in agree- 
ment with the increase (with the laser intensity) of the 
mean number of recollisions required for double ioniza- 
tion. A closer examination of the curve reveals that the 
mean shared action (as well as the mean exchanged en- 
ergy) at recollision decreases as 1/Eq. This decrease is 
also found rigorously using a simplified model of recol- 
lision between an inner electron in a harmonic potential 
recolliding with a free outer electron. 
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FIG. 2: Standard deviation (red stars, left hand y-scale) of 
shared action during recollision for the full Hamiltonian JT]) 
as a function of the laser intensity /. An approximation of 
the standard deviations (red continuous curve) is given by 
Eq. ([3}- The thin black continuous line corresponds to the 
maximum recollision energy £ max [HI, [l3| (expressed in terms 
of actions) and the vertical dotted line indicates the intensity 
after which we stop the mapping. The dashed-dotted curve 
(blue line, right-hand scale) shows the mean number of recol- 
lisions computed from all analyzed trajectories. The typical 
CDI process requires between one (low intensity) and four 
(high intensity) recollisions. 

To combine the two trends of the mean shared action 
AA(Eq) (proportional to Eq at low intensities as given 
by £ max and to I/Eq at higher ones) we fit it by: 



AA(E ) 



aEl 



(3) 



which combines the two features. The parameters a and 
(3 are computed so as to accurately reproduce the evo- 
lution of the mean exchanged energy during recollisions 
(see continuous lines in Fig. here a = 3.5 x 10~ 2 
and /3 = 1.5 x 1CT 4 . 

Mapping model for recollisions. A sequential double 
ionization involves an inner electron which is classically 
confined on the invariant tori of Hamiltonian ^ [19j . 
In the neighborhood of the nucleus, the motion is har- 
monic with a frequency of y/2, and moving away from 
the nucleus, the frequency decreases. We consider a 
simplified model for the inner electron in action-angle 



variables given by the integrable Hamiltonian Hq(A) = 
-2 + V2(e aA -l)/a where a can be adjusted for quantita- 
tive agreement with the frequencies obtained for Hamil- 
tonian 1(3} in the bound region (here, a = —0.8). This 
formula for Ho comes from the observation that the fre- 
quency depends approximately linearly on the energy in 
the whole bound region. 

When the outer electron returns to the core, it gives 
a kick in action to the inner one, which jumps from one 
invariant torus to another (or to the unbound region), 
i.e., at each recollision, the inner electron experiences a 
kick. The action of the inner electron is constant between 
two recollisions. This collision dynamics is modeled by 
the "kicked" rotator (H: 



H m {<p,A,t) = H (A)- 



N 



-eAcosip S (t - nT), 



(4) 



where Hq(A) is the integrable part of the Hamiltonian 
of the inner electron, and T is the delay between two 
recollisions, which is assumed to be constant and equal 
to half a period of the laser field. We denote by N the 
number of recollisions. The recollisions are modeled by 
a kick in action equal to eA cos tp such that a kick might 
increase or decrease the action according to the respective 
phase between the two electrons. In addition, it is more 
difficult to kick the inner electron out if it is at the bottom 
of the well, so the kick strength is proportional to A. In 
this way, the action remains positive at all times. The 
maximum strength of the kick depends strongly on Eq as 
given by Eq. (|3j). 

Let ip n and A n be the angle and action of a trajectory 
of Hamiltonian ((4]) at time (nT)" (right before the n th 
kick). By integrating the trajectories between two kicks, 
we approximate the dynamics of Hamiltonian ([¥]) by the 
two-dimensional symplectic map : 



A n+ i = A n /(1- esinifn), 

fn+l = + (A n+ i) T + £ COS ip n , 



(5) 



where ujq (A) = ^/2e aA is the frequency of the inner elec- 
tron. In Fig. [31 we display two phase portraits of the 
mapping §5§ for two laser intensities : One at low in- 
tensity (I = 2 x 10 14 W ■ cm~ 2 ) in the range of inten- 
sity where CDI is maximum, where the phase portrait 
appears to be very chaotic, and one at high intensity 
(I = 2 X 10 15 W • cm~ 2 ) where the phase portrait is 
more regular. In the chaotic region, the diffusion is much 
stronger at the maximum of CDI than for larger inten- 
sities (see Fig. |31 left panel, where trajectories escape 
quickly from the core region, explaining that there are 
fewer points than in the right panel). Since the strength 
of the kicks decreases with the intensity at high intensi- 
ties, the phase space becomes more regular. If the inner 
electron is inside an elliptic island (which occurs mainly 
at high intensities), it will not ionize regardless of the 
number of recollisions it undergoes. As the intensity in- 
creases, the recollisions become less effective and the map 
becomes integrable so fewer CDI events occur. 
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FIG. 3: Phase space portrait of some trajectories of the map- 
ping (J5]) for low intensity (left panel) I = 2 x 10 14 W • cm -2 , 
and for high intensity (right panel) I — 2 x 10 15 W • cm -2 , 
represented by vertical dashed lines in Fig. [1] In the right 
panel, we indicate the critical action A m = 0.74 by a hori- 
zontal line (whereas the critical action A m is 1.57 in the left 
panel) . 

Through the previous mapping, we have derived a sim- 
ple model for the dynamics of the inner electron initially 
in the bound region which experiences recollisions with 
the outer one. From this model, we compute ionization 
probabilities of the inner electron from which we deduce 
the probability of double ionization as a function of the 
laser intensity: Double ionization occurs if and only if 
the inner electron ionizes and we assume that the outer 
electron remains ionized for all times. The picture of 
the bound and unbound regions for the effective Hamil- 
tonian (|2|) of the inner electron gives a natural criterion 
for ionization. We recall that once the inner electron 
has reached an action larger than the outermost invari- 
ant torus (with action A m ), it is driven away from the 



nucleus by the laser field. Therefore, all recollisions lead- 
ing to an action larger than a critical value A m (which 
depends on E ) subsequently lead to ionization of the in- 
ner electron and, in angle-action variables, the unbound 
region becomes V (E ) = {(p,A) s.t. A > A m (E )}. 

From an initial distribution in angle-action coordinates 
obtained from the microcanonical distribution of inner 
electrons in phase space, we iterate the mapping ([5]) 
a fixed number of times for different intensities (and 
thus different e). In Fig. Q] we display the probabili- 
ties of double ionization as a function of the laser in- 
tensity. We disregard recollisions for intensities larger 
than 2.5 x 10 15 W ■ cm" 2 . This adjustment is motivated 
by the weak probability of recollisions we have detected 
in the data analysis process (see Fig. [2]). We notice that it 
qualitatively reproduces the trends observed in the dou- 
ble ionization yields observed using a statistical analy- 
sis of trajectories of Hamiltonian JT]). In particular, the 
asymmetry in the increase and decrease of the CDI prob- 
ability is worth noting. 

In summary, we have identified key properties of the 
inelastic electron-electron collisions through the analysis 
of recolliding trajectories. We connected our findings to 
a simplified model for the dynamics of the recollisions, 
which amounts to a two-dimensional symplectic map in 
the action-angle coordinates of the inner electron. The 
statistical analysis of the (discrete-time) trajectories of 
this model results in the knee shape for the probability 
of double ionization versus intensity. A proper decompo- 
sition into correlated and uncorrclatcd double ionization 
yields a bell-shape for the correlated process and a mono- 
tonic rise for the uncorrclatcd one. 
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